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Abstract: In our previous work the nonlinearity parameter /jvx, which characterizes non- 
gaussianity in the cosmic microwave background, was estimated for a class of inflationary 
models based on nonlocal field theory. These models include p-adic inflation and generically 
have the remarkable property that slow roll inflation can proceed even with an extremely 
steep potential. Previous calculations found that large nongaussianity is possible; how- 
ever, the technical complications associated with studying perturbations in theories with 
infinitely many derivatives forced us to provide only an order of magnitude estimate for 
Inl- We reconsider the problem of computing f^i in nonlocal inflation models, show- 
ing that a particular choice of field basis and recent progress in cosmological perturbation 
theory makes an exact computation possible. We provide the first quantitatively accurate 
computation of the bispectrum in nonlocal inflation, confirming our previous claim that it 
can be observably large. We show that the shape of the bispectrum in this class of models 
makes it observationally distinguishable from Dirac-Born-Infeld inflation models. 
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1. Introduction 

Nongaussianity in the Cosmic Microwave Background (CMB) is usually characterized by 
the dimensionless nonlinearity parameter, /jvx, which is derived from the bispectrum of 
the curvature perturbation. In the simplest models of inflation f^L is suppressed by 
slow roll parameters and hence one has 1/tvlI *C 1 Nongaussianity can further 

be characterized using the trispectrum, which is also small in the simplest models ||. (See 
also ||-f8|] for further details and see || for a review.) On the other hand, if \/nl\ ~ 5 is 
observed it would rule out the minimal inflationary scenario and favour one of a relatively 
small number of nonminimal constructions. Such models include single field inflation with 



a small sound speed pQ ]. An example is Dirac-Born-Infeld (DBI) inflation [11], whose 
generalization to multiple fields has recently been considered in Another way to get 
observably large Jml from single field models is by having one or more sharp features in 



the inflaton potential [13]. There has also been considerable interest constructing multi- 
field models which can potentially give large nongaussianity using canonical fields and 
smooth potentials [14|. Successful examples include the curvaton model [15, |l(| and hybrid 



inflation |L7| where large nongaussianity may be generated during the preheating phase, 
after inflation [18, O]. (See also [2C] for further studies of nongaussianity from preheating.) 
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The current limit on nongaussianity in the CMB is —9 < Jnl < HI 1 for the WMAP 



five-year data [21] and future missions are expected to be able to probe /nl down to 
order unity [22]. Efforts to construct inflationary models with \/nl\ 3> 1 have recently 
received additional motivation from the claim of a detection of 26.9 < /nl < 146.7 at 95 
C.L. in |23| (see also [24|, which claims a lower but nonvanishing detection, with central 



value Jnl — 50, at the same confidence level). In order to assess the significance of such 
detections it will be crucial to be able to measure Jnl using independent techniques, such 
as through its effect of the abudance and clustering of virialized objects In the case 
that /iVL has strong scale dependence cluster number counts may, in fact, provide better 



propects for detection than CMB observations [26]. 



In (27]] a class of inflationary models was constructed based on nonlocal field theories 
containing infinitely many derivatives (see also ^8|). Thep-adic inflation model ]^9| belongs 
to this class. In [^7]] it was shown that nonlocal inflation models can give rise to observably 
large nongaussianty and, for p-adic inflation, it was estimated that /tvl ~ 120 (for a string 
coupling of order unity). However, the calculation of pTl] only considered the perturbations 
of the matter field, rather than also including metric perturbations self-consistently. This 
approximation was necessary in order to circumvent the technical complications associated 
with performing cosmological perturbation theory using Lagrangians with infinitely many 



derivatives. This neglect meant that the calculation of [27] only provided a rough estimate 
of the size of /jvx and gave no information about the shape or running. Furthermore, the 
presence of an infra-red (IR) divergence in the momentum integrals lead to concerns that 
the curvature perturbation might not be freezing in on large scales, as would be expected 



in models with only adiabatic perturbations [30] 



In this paper we reconsider the calculation of the bispectrum in nonlocal inflationary 
models, making use of recent advances both in cosmological pertubation theory [31] and also 



in the formal mathematics of infinite order partial differential equations [^] . The analysis 
is rendered tractable by using a particular choice of field basis and by taking advantage 
of the hierarchy e <C \r]\ among the slow-roll parameters in hill-top inflation models. The 
calculation is further simplified by use of the recent formalism for computing bispectra 



directly from the field equations developed by Seery, Malik and Lyth | 3lJ| . Our analysis 



confirms the qualitative estimate of |27|] while improving significantly on the accuracy. 

The organization of this paper is as follows. In section || we discuss the class of nonlocal 
theories under discussion and demonstrate how the Lagrangian may be simplified by use 
of appropriate field redefinitions. In section || we write out the perturbed Klein-Gordon 
equation at first and second order in perturbation theory. In section |] we use the results of 
[31] to construct the nonlinearity parameter for nonlocal inflation. In section |5| we specialize 



to the case of p-adic inflation, confirming previous claims that /tvl can be made large in 
that model. Finally, in section || we conclude. 



This is assuming the "local" form of nongaussianity. For the equilateral model one has the slightly 
weaker constraint — 151 < /jvl < 253. 
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2. Nonlocal Hill- Top Inflation 



2.1 Lagrangian and Field redefinitions 

We are interested in a class of nonlocal scalar field theories of the form 



7 



(2.1) 



where, following [27], we take 7, m s to have dimensions of energy while 4>, U(<f>) and F(z) are 
dimensionless. Lagrangians of the form ( |2.1| ) are ubiquitous in string field theory |[3| and 
other stringy constructions [34]-[35| and can also be motivated in certain braneworld models 
1 36] or as a low energy description of unparticles [37|. Phenomenological constructions of 
the form (2.1) have also attracted interest due to possible implications for the cosmological 
constant problem y3q] . (Similar nonlocal constructions are expected to arise generically in 
quantum field theories with a minimal length scale 



].) In (2.1) the mass scale m s gives 
the scale of the underlying nonlocality of the theory. For models derived from string theory 
this will coincide with the string scale x/c7, however, we do not require this assumption in 
our analysis. 

We assume a hill-top [4C] form for the potential 



U(<f>) 



TT ^ A? 1 3 ,3 , 

~ 21^ 3T 



(2.2) 



so that (j) = is an unstable maximum. In (2.2) • • • denotes quartic and higher order 
terms which stabilize the potential at large values of eft. Since we are primarily interested 
in dynamics near the unstable maximum we can neglect the terms 0(</> 4 ) and higher. 

In order for the theory to be free of ghosts and Ostrogradski instabilities we require 
that F(z) can be cast in the form 



7 



F I ] + //" 

mi 



□ 



(2.3) 



where T(z) is an analytic function of the complex variable z having no zeros at finite z 
1 32 ] . Notice that, as in [27], we have 



= -F(-u 2 /m 2 s ) 



(2.4) 



and F(0) = 0. We will be particularly interested in the following special values of the 
function T(z) 



r 



= —F' 

7V 2 



r(o) 



uJ- 



The structure (|2.3|) for the kinetic function F 



(2.5) 
(2.6) 

is required in order to ensure that 



the propagator has only a single pole and the resulting quantum theory (2.1) describes 
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only a single degree of freedom. This constraint also ensures that the Cauchy problem 
associated with the classical equation of motion is well-posed with only two initial data 
1 32 J. (See also [41] for an alternative approach to counting initial data.) Note also that 
theories of the form (^^) which do suffer from the presence of ghosts (or classical higher 
derivative instabilities) may nevertheless be rendered physically sensible by appropriately 
redefining the action pseudo-differential operator F ( -3? ) using the prescription advocated 



in ||. 

The equation of motion derived from (|2.1| ) is 



□ 



m: 



u'{<t>) 



where Ui n t denotes the terms in U(4>) which are of higher order than quadratic 

9 a3 



U- u 



3! 



+ 



(2.7) 



(21 



In order to study perturbations it will be convenient to introduce a "dressed" field (p 
(see p2|) defined through a nonlocal field redefinition 



if = 7 r 



□ 



w- 



1/2 



(2.9) 



We will also refer to ip as the canonical inflaton for reasons which will become obvious 
shortly. Such a nonlocal field redefinition is allowed since T(z) contains neither zeroes nor 
poles hence there is no danger of altering the structure of the initial value problem 
]32[ |. In terms of the dressed field the Lagrangian ( |2.1| ) takes the form 



£= {a + uj 2 ) if - V - V int 



-1/2 



where Vq = 7 4 C/q and 



Mnt [X] 



(2.10) 



(2.11) 



When □(■ 



-to 



(as occurs in the linear theory) then ( |2.9D reproduces the definition of a 



"canonical" field employed in |28| although it differs from the definition used in [27]. The 
equation of motion derived from ( p. 10 ) takes the form 



(□ + u; 2 )ip = r 



□ 



m- 



-1/2 



□ 



-1/2 



(2.12) 



2.2 Background Evolution 



In [£7J the homogeneous solutions of the theory (|2.l| ) were constructed using an expansion 
in powers of e xt . In [29] the backround solutions of p-adic inflation were constructed using a 



similar expansion and were verified using an independent formalism (the friction-dominated 
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approximation) which did not rely on small values of 0. The result is that the background 
canonical field <po(t) and the Hubble scale H(t) evolve as 



where 



<p (t) Ae^ Hot 

H(t) ^ H - H 2 e 2 ^ Hot 



Vo _7 4 ^o 
Ml Ml 



7 



mi 



(2.13) 
(2.14) 



(2.15) 
(2.16) 



(using eq. (E 



and eq. ( |2.5| )) and in (2.13- 2~l4| ) we have defined the slow roll parameter 

,2 



f] 



UJ 

ml 



(2.17) 



(The explicit expression for H 2 is given in [27|, H 2 = uj 2 HqF' y—^J /[4Uqiti 2 ].) Inflation 
takes place at large negative values of t so that e^ Hot <C 1 throughout. The zero of time 



is chosen abitrarily in the solution (2.13), but this of course has no effect on any physical 
observable. We define a second slow-roll parameter 



H 
W 



2MIH 2 



(2.18) 



The remarkable feature of the solutions ( |2.13| , 2.14 ) is that slow-roll inflation proceeds even 
when the naive potential (which would be inferred by setting □ = in (|2.1[ )) is extremely 
steep. 

Notice that, as is typical of hill-top models, we have large hierarchy between the slow 
roll parameters e <C \r)\. It is straightforward to verify that 

\m h 

so that in the limit e/\r]\ — > one has e 2 ^ Hot -* and thus <p (t) -> 0, H(t) -* H . 
Following [^] we conservatively take the end of inflation to be when e ~ |r/|, rather than 
e ~ 1, since this is the point where the perturbative expasion ( 2.14| ) breaks down. 2 (In the 
case of p-adic inflation the friction dominated approach gives the same result for the end 
of inflation p9[ .) It follows that it is consistent to work with a pure de Sitter background 
<Po(t) =0, H(t) = Hq while keeping terms up to leading order in the rj slow-roll parameter 
to study perturbations. 



2 The 0(e 2MHot ) correction to fl2.15| ) was computed in 
subdominant to the leading term even when e > |?j|. 



where it was shown that this term is typically 



- 5 - 



3. The Perturbed Field Equations 



3.1 Linear Perturbations 

The computation of the bispectrum for the theory fl2.1|) is greatly simplified by using 
the formalism of Seery et al. for working directly with the perturbed Klein-Gordon 
equation. Since ( 2.1C| ) has canonical structure in the limit V^t — ^ it is clear that ( [2.12 ) is 
the appropriate starting point to match onto the analysis of |3l]]. We perturb the field as 



(p(t,x) = (po(t)+S(p(t,x) (3.1) 



and work in uniform curvature gauge (neglecting tensor modes) so that the ADM [43] 
metric can be written as 

ds 2 = -Ndt 2 + 5ij (dx i + N^t) (dx j + N j dt) (3.2) 

It will be convenient to work in terms of conformal time r defined by adr = dt so that, 
in the limit e — ► 0, one has a = —1/(Hqt). We denote derivatives with respect to r by 
/' = d T f and introduce a conformal time Hubble scale 7i = a' /a = a. 
Following we decompose the field perturbation as 

6(p(t, x) = 8i(p(t, x) + -6 2 cp(t, x) (3.3) 

where 5np is defined so that it obeys exactly gaussian statistics. Since the nonlocal structure 
of Vint i n the Lagrangian (2.10) does not appear at the linear level, the equation for 8\tp is 
precisely the same result that one would have in a standard (local) model 3 , derived in [45] 



a 2 V 

Sw" + ZHSitft - V 2 5np - a 2 J 2 5^ - 2e — T 5 1 ip = (3.4) 



The final term in ( p.4[ ) arises from the terms in Hip involving the metric perturbations 
which, to leading order in perturbation theory, always involve time derivatives of (fo and 
hence are proportional to e. In the limit e — > these terms vanish and ( p.4[ ) is identical to 
the dynamical equation for a light field in de Sitter space. Denoting the d'Alembertian of 
de Sitter space by Do, that is 

□ = -d 2 t - 3H d t + e - 2Hot V 2 (3.5) 
we see that fl3.4|) is equivalent to 

□ <5 1(/ 9 ^ -uo 2 bxip (3.6) 

for e — > 0. In other words, in the limit e <C we can safely replace n5\ip with □o^iy 
throughout the calculation. This trivial observation will greatly simplify our treatment of 
the nonlocal structure of the theory ( 2.1 Q| ) . 



3 In |$^] the equivalence of (^^) to a local theory at the linearized level was demonstrated using the 
formal calculus of pseudo-differential operators. 
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The solutions of (|3.4j ) are well-known 



d 3 k 



ikx 



(2vr) 3 / 2 

6k(t) = «k^k(i) + at. k p*_ k (t) 



7T 



2 V a 3 tf v \aH 
where the order of the Hankel function is 



Ll) 



4 + m^2 + OM 



(3.7) 
(3.8) 
(3.9) 



and we have normalized the annihilation/creation operators as 
Q3.9j ) the constant real-valued phase 5 is irrelevant for our calcu 
the rj parameter we have 



e- ikT {\ + ikr) 



flk, 4, =5®(k-k'). In 
ation. To leading order in 

(3.10) 



Using the solution ( |3.7|]3.9| ) it is straightforward to construct the curvature perturba- 
tion in the linear theory. The analysis is essentially identical to [27], the only difference 
being that ( |2.9| ) implies a slightly different normalization for the canonical scalar which, 
in turn, gives a slightly different expression for the COBE normalization. In the appendix 



we redo the linear analysis of [27] taking into account this change of normalization. 



3.2 Second Order Perturbations 

We now turn to the equation for 8 2 (p. All of the contributions coming from the left-hand- 
side of (2.12) appear already in the standard (local) theory and have been computed in 
[44, 45 1 . The only new contribution coming from the nonlocal structure of ( [2.10D arises 



due to the term V\ n t [T 1 / 2 (U/rn 2 s )Lp\ and appears on the right-hand-side of ( 2 . 1 2| ) . Thus, 
at leading order in slow roll parameters, the equation for 8 2 ip is 



6 2 <p" + 2H5 2 ^ - V 2 5 2 <p - a 2 uj 2 5 2 cp = J in t,2 + F 2 (S 1 <p) + G 2 (5np) 



(3.11) 



where the source terms F 2 , G 2 are the same as those appearing in [31]. The quantities F 2 , 
G 2 were explicitly computed in terms of field perturbations in [44, ^] and their contribu- 
tions to /jvx were derived in [j3lj. The only new term in (3.11) is J\ n t,2 which denotes the 
leading contribution in a perturbative expansion of 



Ji 



mi 



mi 



-1/2 



□ 



???; 



-1/2 



(3.12) 



To evaluate J- m t. 2 we first consider the argument of V^ nt in ( 3.12| ). As we discussed 
above, in the limit e <C \r]\ we can take cpo = and safely replace U5\(p by O 6np (where 
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□o is the d'Alembertian in de Sitter space, given by (|3.5|) ). Thus we have 



2 \ -V 2 

' 5iLp (3.13) 



m 2 



where the corrections are either proportional to e or are higher order in perturbation theory. 
On the second line of ( |3.13j ) we have used the fact that □o^iV? — ~~ uj 2 5\(p in the small-e 
limit (see eq. pT6| ). 

Using ( 3 . 1 3| ) and focusing on the cubic term V\ n t(x) = g^/x 3 /6 we find that 



in; 



J^V^r ffl (5 W ) 2 (3.14) 



/ \-l/2 

The challenge is now to compute the action of the nonlocal operator r ( ^-J on the 
product (5iip) 2 . This is, in general, extremely complicated since (^i^) 2 is not an eigen- 
function of Dq. However, since the perturbation 5 up behaves as a light field in de Sitter 
space, the result will be fairly simple. To see this, we explicitly compute the expectation 

value (r(Vf 1/2 (<M 2 >. 



For any analytic pseudo-differential operator g(O ) we have 

___ e *k.x 9M 2 _ 3Hodt _ k 2 /a 2 } {&m l_ k ( t )) 

^fe^g(-df - 3H d t - k 2 /a 2 ) 6® (k) \p k ,(t)\ 2 

72' 



(2vr) 3 2/fc 3 V a 2 H 2 g(0) 



(3.15) 



Thus, we have the identity 

(<7(n )(M 2 (*)> = 5(0) ((5np) 2 (t - At)) (3.16) 

where the time shift is 



1 , 
At = In 

2H 



9(2H 2 ) 



(3.17) 



5(0) 

In ( 3.16j ) the overall factor of g(0) arises due to the large-scale constancy of S\cp while the 
time shift At reflects the fundamental nonlocality of the underlying operation. Although, 
in general, such a time shift could have a significant impact on the calculation it turns out 
to be unimportant here because the variance ((5ip) 2 } ~ Hq is very close to a constant. 
Thus, we can safely approximate g^oX^i^) 2 ~ g(0)(8np) 2 . This result is intuitively 



-8- 



reasonable because (5iip) ls approximately constant on large scales while on small scales 



it is negligible. Setting g(O ) 



□o 



-1/2 



we have a time shift 



At 



4fln 



In 



r(o) 



(3.18) 



which amounts to less than a single e-folding in the case of p-adic inflation. 
In summary, we find that 



91 



int 



r(o)V2 



(3.19) 



which takes exactly the form which would arise from a cubic interaction of the form Vint 
cy(<5i<£>) 3 where the effective cubic coupling is 



Cy = + 



91 



1 



3! 



r(o)V2 



(3.20) 



Equations and (^20|) are the main results of this section. 

Notice that ( |3.20[ ) does not exactly coincide with the definition of the quantity cy in 
|27j, due to the more quantitative treatment we are carrying out here. However, for the 
interesting case of p-adic inflation both results are the same order of magnitude. 



4. The Bispectrum and Nongaussianity 
4.1 Results for f^L 

The result ( |3.19| ) is extremely simple. Since a term of this form was already included in [31] 
it follows that we can obtain /tvl simply by substituting V' = 0, V" = —uj 2 and V'" = 6cy 



into the results of [31]. As in [27], large nongaussianity will result due to the fact that the 
effective cubic interaction can be made extremely large without spoiling slow roll. This is 
a result of the nonlocal structure of the theory; the higher derivative structure in V\ n t can 
mimic a strongly interacting theory even when the actual couplings are all perturbative. 
(From the perspective of the undressed Lagrangian (|2.1| ) one would say that the nonlocal 
corrections to the kinetic term conspire to slow the rolling of the inflaton even when U((f>) 
is extremely steep.) Notice that the size of cy is, however, bounded from above by the 
requirement of perturbative stability against quantum corrections. This bound is discussed 
in 0. 4 

In order to make contact with the analysis of |3l]] we must construct and effective 
value for the slow-roll parameter £( 2 ) = Mf V'V'"/V 2 which parameterizes the size of V" 
in a local inflationary model. In our case, however, £( 2 ) can greatly exceed unity without 



violating the conditions for slow roll; we will discuss this issue in section 4.2 . 



4 For the case of p-adic inflation the bounds of 
one obtains in the full string theory set-up. 



are equivalent to g s 1 which is precisely the result 
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(2) 

To compute an effective value of Q ff in the nonlocal model, which will play the same 
role relative to /jvx as does £( 2 ) in local models, we need to replace V" by an effective 
value V^g, which can be deduced by comparing our result for Ji nt , eq. (|3.19| ), with ref. |H]]'s 
corresponding expression, J; nt = —a 2 V"'5i(j) 2 . This gives = 6cy. To compute V in 
conventions of ref. |3l|], we should work in terms of the canonically normalized field (p, so 
that V = —ui 2 ip = — u 2 ^q/ {\rj\Ho) = — (PqVq/{M 2 Hq). In this way we find that 

V" 

- -sign^v^M 3 -^ (4.1) 
where the slow roll parameter e, defined in ( 2.18| ), is understood to be evaluated at horizon 



crossing. Equation (|4.1|) can be rewritten in a form similar also to that given in ref. \21\ 



where is defined by (|A-7|) and we have used ( A-9[) . (Equation (|4.2|) should be compared 
to eqn. 4.15 in 



(2) 

In order to give a sense of how the quantity scales with the various couplings in 



the theory we substitute the explicit expression for cy (13.201) into (fO) and use (12 . 5|-pT3) 

— — — 

as well as the results of the appendix to rewrite Q s in the form 

Cff - -2ffsign(<^o)— tt^ (4.3) 

(2) 

This expression shows explicitly how Qs depends on the model parameters /i, Uq, g, m s 



F(z) and the derived "effective mass" u. However, eq. (4.3) does not take into account 
the various constraints among these quanties which are required for consistency with CMB 
observations. 

(2) 

To more explicitly evaluate P ff , we would like to incorporate the experimental con- 
traints on the amplitude and index of the CMB power spectrum. These constraints allow 
us to fix two combinations of the model parameters, which can be chosen as (i 2 and 7 4 C/o- 
It is convenient to express these in terms of the combination 

f A = 127t 2 Ale- N ^ n °- ll \n s -l\ 2 ~ (0.1 - 1) x 1(T 10 (4.4) 

where A 2 = 25 x 10~ 10 is the COBE normalization for the amplitude of the scalar power 
spectrum (see the appendix for more details). The quantity Ja in ( |4.4j ) depends upon 
observables rather than parameters of the model. In figure |l| we show the dependence of 
Ja on n s for a range of iV e . Using eqs. ( |2.4j ), QA-5 ) and ( |A-9| ) , we find 




f = ~F \~fA-^ )^^ = M$f A (4.5) 
7 4 [/ = 2M p 4 /aK " (4-6) 
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0.9 0.92 0.94 0.96 0.98 1 

n s 



Figure 1: /a versus n s , for N e — 40, . . . , 60. 



Using these equations to eliminate \j? and Uq, we find that 

-sign(</?o) 



t(2) 
Seflf 
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f A M p r{-f A M 2 /m 2 ) V2r(0)e^l" 



1/2 



(4.7) 



(2) 



The expression for the nonlinearity parameter can be written simply in terms of 



as |i| 



NL 



e(2) 
s e ff 



AT e + 0.91 + [ k t hkj - -fcf 



i<j 



9 



+ ... 



(4.8) 



where kt = ki and • • • denotes the (subleading) terms which appear in the standard cal- 
culation. As was explained in [31], the dependence of the number of e-foldings since horizon 
crossing A^ does not imply super-horizon evolution of £ since this explicit dependence is 
cancelled by implicit N e dependence in the parameters and r/. 

To estimate the size of the nongaussianity we evaluate (|4.8| ) on equilateral triangles 
ki=k 2 = k 3 



J NL 



6 



(A r e - 2.09)C, 



(2) 
off 



(4.9) 



Equations (4/7|- 0|) are the main result of this section. Using these expressions one may 
straightforwardly evaluate Jnl in terms of the model parameters g, Uq, [/,, m s , subject to 
the constraints (4J3-0). Note that w is a derived parameter, u 2 = — m^i ?_1 (— fi 2 ). We 
will confirm that, unlike in standard single-field slow roll inflation, /jvx can be quite large 
in nonlocal models. 

Finally, let us comment on our choice of sign convention for Jnl- Our definition 
follows and coincides with the convention used by the WMAP team [21| and also by 



5 One may verify using the results of [^] that the relevant flow equations continue to hold even when 
is larger than unity. It follows that the same arguments given in |u| imply that £ is frozen on super-horizon 
scales in nonlocal inflation, at least to the accuracy of our calculation. 
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1 23], however, it differs from some other studies || [lC]]. Notice that the final result for f^L 
in []3l| actually has a sign error for the terms involving e* and 7/*, however, this error does 
not affect the £^ term coming from the cubic interaction. 6 Since it is this latter term that 
dominates our result fl4.8| ) it follows that our results can be directly compared to pi] , g3fl . 
For a careful discussion of sign conventions for /jyx the reader is referred to appendix A. 2 
of Pi. 



4.2 Consistency of the Slow-Roll Dynamics 

Even in a conventional local model, one could always arrange for V" to be anomalously 
large at some special field value (p c , and seemingly obtain a sizeable contribution to the 
nongaussianity. However this would not be consistent, because the conditions for slow roll 
would be spoiled for <p ^ <p c , leading to a rapid end of inflation, long before the required 
number of e- foldings. We need to verify that the same problem does not occur in the 
nonlocal models under consideration. This is a somewhat subtle issue, and one can get the 
wrong answer by being too naive. 

(2) 

The crucial question which we must address is whether the large values of required 



to get large /tvl will invalidate the slowly-rolling backgound solutions (|2. 13 -2.14). In the 



case of p-adic inflation we have verified both the background solutions (|2.13 - 244]) and also 



the expression for the end of inflation e ~ \rj\ using the friction-dominated approximation 
|29f| . This formalism is nonperturbative and gives a powerful consistency check on the slow 
roll dynamics of p-adic inflation in the |/jvx| ^ 1 regime. 

However, for a general model of the form fl2.1|) it may be nontrivial to verify the 
consistency of the small e'* 7 ' expansion in d2l3tf03 ) when \f NL \ > 1. This is so 



because |/jvl| 3* 1 tends to require \g\ 3> 1 and it is conceivable that terms involving 
higher powers of e^ H °\ w hi c h were neg lected in ( glfjng ), become important before 



e ~ \rj\. We caution the reader that, although we present our results in the most general 
context, the consistency of the slow roll dynamics should be verified for any particular 
model under consideration, as we have done for p-adic inflation. 



5. p-adic Inflation 

The Lagrangian of p-adic string theory is 



9 

9p 



■n/2m 2 s 



p + 1 



ll'- 



p+1 



(5.1) 



with gp = g1{p — l)/p 2 and one must take ij) = 1 + (j> to put this in the same form as (|2.1|). 
The constants g s , m s are the string coupling and string mass respectively and p is a prime 
number that characterizes the field of numbers in which the worldsheet coordinates take 
values. (A p-adic number q can be expressed in the form of a series q = YlnLk c nP n where 



6 We are grateful to D. Seery for clarification on this point. 
7 Note that ref. [^7| missed the factor of 2 in the exponent p~ D//2 " 
propagated into the final results. 



This typographical error was not 
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k is a possibly negative integer, k > — oo.) In [27] expressions for the parameters Uq, g, to, 



fj,, ■ ■ ■ in terms of p, g s and m s were derived. They are given by 

M 2 =P"1, ^ = 2ml l A = j^- 2 , U = 1 p±- g = -\p(p-l) (5.2) 

(P - l )9s 2 {P + 1) 2 

and the kinetic function is .F(z) = 1 — p~ z / 2 . Using the results (2.5,2.6) we obtain the 
relevant values 

w 2/ 2^ P ln P P 2ln P w n \ PVP- 1 / KQ \ 

r(-w /mj = — — = 7==, r(0) = — (5.3) 

2c/ p 2g s ^p-l 2g s 



The experimental constraints ([4.5-4.6) then show that 



f 



rl/2 

V2^ P ' 
(P+1)5 S 2 1 



p 2 16 



77t/aK - 1| (5.4) 



We will use (5.4) to eliminate m s and g s , leaving p as the only free parameter. Using the 



values (|5.3|) in (J4.T|) we find that (again for equilateral triangles in fc-space) 

fL = 5(A V 2) K - 1| V*^ (5.5) 
JNL My/2 ' ' lnpv^+T V ' 

which remarkably has no dependence on the small parameter /a, and which scales as /jyx ~ 
p 1 / 2 /^^ in the interesting region of parameter space p> 1 where the nonlocal structure 
of the theory is playing an important role in the dynamics. (This is slightly suppressed 



compared to the result of [27| by the factor hip although the order unity prefactors are 
slightly larger so that the net result has a very similar magnitude.) In the limit p — ► 1 the p- 
adic theory (|5.1|) reduces to a local field theory and (|5.5|) gives |/jvx| ~ \n s — 1|, as expected. 
But in the opposite regime of p 3> 1, we can obtain large nongaussianity, as illustrated in 
figure ^. We have also evaluated Jnl as a function of g s for the WMAP5-preferred value 
of n s = 0.96, and the range N e = 40 — 60 e-foldings, plotting the result in figure |3|. The 
curves are approximately linear, and can be well-fit in the region of observational interest 
by 

/^ L ^(23 + 6.767V e )<7 s (5.6) 

It shows that for reasonable values of the string coupling, 0.15 < g s < 0.25, we reproduce 
the central value for Jnl claimed by ref. [23|. The corresponding values of p ~ 10 12 as a 



function of g s from eq. ( |5.4| ) are shown in figure |4|. The large value of p (as long as g s is 
not extremely small) is a consequence of the smallness of /a ~ A 2 which enters through 
the COBE normalization in (|5.4j), A similar situation occurs in the DBI inflation models 
1 11], where the dimensionless 't Hooft coupling must be taken to be A ~ 10 12 <? s . 

It is worth noting that the large /nl in p-adic inflation is not coming solely from 



the factor N e in (|4.8|) , Indeed, for g s ~ 1 we still have Jnl ~ 8 - almost four orders 
of magnitude larger than in a local field theory - taking only N e = 3. Rather, the large 
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400 



< Z200 




Figure 2: f NL for equilateral triangles, as a function of n s , for the same values of N e = 40, . . .60 
as in figure [j], and values of the string coupling g s — 0.05, 0.15, 0.5. 

250F ^ ^ ^ 




Figure 3: f NL versus as a function of g s , for n s = 0.96 and N e — 40, . . . 60. Horizontal lines show 
central value and 95% c.l. range for detection claimed by ref. B3|. 



nongaussianity arises purely due to the fact that the potential U (eft) is quite steep and the 
cubic coupling is large when p is large: |<7| ~ p 2 3> 1 in the p-adic model. 

An interesting question is whether the result (4.8) can be observationally distinguished 
from the prediction for models such as DBI inflation. In ref. ]10| it was shown that the 
DBI model predicts a shape proportional to 



(5.7) 



where S = Y^i > which differs from the shape predicted for the nonlocal models, (|4.8| ). Let 
us call the latter Anl, and normalize A c and A^l such that they are equal for equilateral 
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Figure 4: p versus g s for n s = 0.96 and N e = 40, . . . 60. 



triangles [k\ = ki = ks). A measure for how easily the two shape dependences could be 
distinguished from each other is the quantity 

Both A c and Anl depend only on the shape and not the size of the triangles defined by k\, 
ki and ^3 = — (ki + ki)- Thus we can take R to be a function of any two angles, call them 
ai 2< R also depends on N e via A^l- Setting fe 3 = 1, we can write ki = sinaj/ sin(ai + 0-2). 
The physical values of oti are within the triangle shown in figure |5[ whose boundaries 
represent degenerate limits: if either ai vanishes, then two of the vertices coincide, while 
if a.\ + «2 = then all three vertices lie on a line, but do not in general coincide. We plot 
contours of R in figure [5] for N e = 55, but we find that their shape is quite insensitive to the 
value of N e , in the range 40 to 60, due to the fact that N e dominates over the /c-dependent 



terms in eq. (4.8), so that Anl is nearly independent of angles. On the other hand, A c 
shows strong dependence on the angles, and the fact that R deviates significantly from 
zero over a large range of angles indicates that the search for this shape-dependence could 
be a strong discriminator between nongaussianity coming from nonlocal inflation and that 
coming small sound speed models inflation models. 

Using the bispectrum alone it may, however, be more challenging to discriminate obser- 
vationally between nonlocal inflation and the curvaton model p5|, since the latter predicts 
/nl of the same form as the leading term in Q4,8| ) (proportional to _/V e ) A similar 

analysis to the one described above for DBI inflation shows that the quantity R, which 
disciminates between nonlocal inflation and the curvaton model, will be at the percent level 
when Ae ~ 30 — 60 (R could be larger for smaller values of N e ). However, we note that if 
|/tvl| 3> 1 is observed then it may be much easier to distinguish between the two models 



using other constraints, such as those coming from isocurvature modes [48] which should 



be significant in the curvaton model when \/nl\ ^ 1 but are absent in nonlocal inflation. 
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Notice that, because we have normalized Anl = A c on equilateral triangles, it is 
irrelevant for our comparison that uses a different sign convention for /jvl than we 
do: the quantity R only contains meaningful information about relative shapes of the two 
bispectra but not the amplitude or sign. 



6. Conclusions 

In this paper we have reconsidered the calculation of the bispectrum in nonlocal hill-top 



inflationary models. We have confirmed previous claims [27] that such constructions can 
give rise to observably large nongaussianity in the CMB while significantly improving on 
the accuracy of previous estimates. Large nongaussianity seems to be a generic prediction 
of nonlocal hill-top inflation in the sense that /jvx tends to be large whenever the nonlocal 
structure of the theory is playing a significant role in the inflationary dynamics. 

We have shown that nonlocal hill-top inflation gives a distinctive prediction for the 
shape of the nonlinearity parameter which should be observationally distinguishable from 
other models which give |/jvx| S> 1, such as DBI inflation (or similar models where large 
nongaussianity is generated due to a small sound speed for the inflaton) . It is more difficult 
to observationally distinguish between nonlocal inflation and the curvaton model using only 
the bispectrum; however, it might be possible using constraints on isocurvature modes. 
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In the interesting case of p-adic inflation |29| we have shown that the predicted value 
of Jnl can be quite close to the central value of the claimed detection in p3[| , using natural 
values for the string coupling. Indeed, future missions should be able to robustly rule out 
(or confirm!) the most interesting region of parameter space for this model. Since p-adic 
string theory is not construed as a realistic model by itself, it remains an open problem to 
construct similar nonlocal inflationary models in more realistic string theory constructions 
(addressing also issues of compactification and moduli stabilization). On this front we are 
optimistic because the nonlocal structure (|2.1|) is ubiquitous in string field theory and can 
be motivated also in a variety of other contexts. 

From the formal perspective, one of the most interesting features of theories of the 



form (2A) is that the both the consistency of the underlying quantum field theory, and 
also the prediction \/nl\ S> 1 in the inflationary context, rely on UV completion. In other 
words, if the Taylor series expansion of the kinetic function F{z) is truncated at some 



finite order 0(z N ) then the resulting truncated theory will, in general, contain spurious 



ghost-like degrees of freedom which are not present in the full theory (2.1) defined by the 
exact transcendental expression for F(z) p2| . Moreover, the dynamics of the truncated 
theory may be quite different from those of the fully nonlocal theory and, in general, our 
prediction for Jnl j ) will not be reproduced. Remarkably, this statement continues to be 
true even for N 1. One might consider the requirement of UV completion daunting were 
it not for the fact that constructions of the form ( |2.1j ) can be microscopically justified from 
string theory. From this perspective we have the exciting possibility that observations of 
the primordial bispectrum might be used to probe physics in the quantum gravity regime. 
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APPENDIX: Constructing the Linear Curvature Perturbation 



In this appendix we redo the linear computation of the curvature perturbation in 27] taking 



into account the change of normalization implied by (2J3). Since, for the linear equation 
( |3.6| ) we have □o^iv? = — uj 2 5\ip equation (|2.9|) simply gives 



6i<p = A6 1( j> (A-l) 
where, using the identity (|2.5[), we have 



,2 _ „, 2t , / U 2 \ _ OJ 2 



^=7 2 r(-^)=^F'(-^) (A-2 



from ( |2.9| ) and (pT5|). The power spectrum associated with the solution (3.7-Q|) is well- 
known 



p *=^ bar' (A - 3) 
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with spectral index 



1 = 2i] 



2uj z 
"3flg 



so that n s < 1, in agreement with the latest WMAP data 
3Hq = 7 4 C/o/Mp equation (A-4) implies the useful relation 



1| Co 



7* 



M p 2 u; 2 



(A-4) 

Using the fact that 
(A-5) 



We now construct the curvature perturbation. In the comoving gauge we have [44] 

C = -—5 W (A-6) 

Following pi we define the constant 

c c = -- (A-7) 



with the understanding that the right-hand-side is evaluated at horizon crossing. The 
COBE normalization fixes the amplitude of the power spectrum of the curvature per- 
turbation 



C 



A 2 



k 



aH 



n a — l 



to be A 2 . = 25 x 1CT 10 . From flA^gp and ( gg ) it is clear that we have the relation 



We now evaluate ( |A-7| ). Using the background solution of j27]] we have 



H 2 1 



^ c A 2 \ 2 ul 



(A- 



(A-9) 



(A-10) 



where u = e xt , A = \r]\H$ and denotes u{t) evaluated at the time of horizon crossing. 
Noting that u(t) = a(t)^ and \r)\ = \n s - l|/2 we have 



= e 2 



■K-i| 



(A-ll) 



where N e is the number of e-foldings between horizon crossing and the end of inflation 
(typically N e ~ 60). We conservatively take inflation to end when e ~ \r]\ since this is 
when the expansion in powers of u(t) becomes invalidated. This gives (eqn. 2.28 of |27j) 



— I F' ( -— 7T 



uJ 



m- 



(A-12) 



Plugging flA^), flA4j) , (|A-11|) and (|A-12[) into (|A-10|) we derive the useful identiy 

° c ~ 7 4 C/ K - 1| 2 



(A-13) 
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Using ( A-13 ) and the COBE normalization ( A-9 ) we get an expression for m s /M p 



'"- ^ = 121r 2 A 2 c \n x -l\ 2 e- N ^ 1 ^ 



M, 



(A-14) 



For the particular case of p-adic inflation equations ( |A-4| ) , (|A-9| ) and ( [A- 10 ) imply the 
following relation between p and g s 



4 gg(p+ 1) e^l""" 1 ! 
37r 2 p 2 \n s — 1| 3 



C 



(A-15) 



which is almost identical to eq. (5.20) in [29|. It follows that this change in normalization 
is not quantitatively significant in the case of p-adic inflation. 
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